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ABSTRACT 

Satellite galaxies in rich clusters are subject to numerous physical processes that can significantly 
influence their evolution. However, the typical L* satellite galaxy resides in much lower mass galaxy 
groups, where the processes capable of altering their evolution are generally weaker and have had less 
time to operate. To investigate the extent to which satellite and central galaxy evolution differs, we 
separately model the stellar mass - halo mass (M* — Mh) relation for these two populations over the 
rcdshift interval < z < 1. This relation for central galaxies is constrained by the galaxy stellar 
mass function while the relation for satellite galaxies is constrained against recent measurements of 
the galaxy two-point correlation function (2PCF). Our approach does not rely on the abundance 
matching technique but instead adopts a flexible functional form for the relation between satellite 
galaxy stellar mass and subhalo mass, where subhalo mass is considered at the maximum mass that a 
subhalo has ever reached in its merger history, M pca k- At z ~ the satellites, on average, have ~ 10% 
larger stellar masses at fixed M pca k compared to central galaxies of the same halo mass (although the 
two relations are consistent at 2 — 3er for M pca k > lO 13 Af ). This is required in order to reproduce 
the observed stellar mass-dependent 2PCF and satellite fractions. At low masses our model slightly 
under-predicts the correlation function at ~ 1 Mpc scales. At z ~ 1 the satellite and central galaxy 
M* — Mh relations are consistent within the errors, and the model provides an excellent fit to the 
clustering data. At present, the errors on the clustering data at z ~ 2 are too large to constrain the 
satellite model. However, we find that a simple model in which satellite and central galaxies share 
the same M* — Mh relation is able to reproduce the extant z ~ 2 clustering data. We speculate that 
the striking similarity between the satellite and central galaxy M* — Mh relations since z ~ 2 arises 
because the central galaxy relation evolves very weakly with time and because the stellar mass of the 
typical satellite galaxy has not changed significantly since it was accreted. The reason for this last 
point is not yet entirely clear, but it is likely related to the fact that the typical ~ L* satellite galaxy 
resides in a poor group where transformation processes are weak and lifetimes are short. 

Subject headings: cosmology: theory — dark matter — galaxies: halos — galaxies: evolution - 
galaxies: clustering — large-scale structure of universe 



1. INTRODUCTION 

Understanding the evolution of stellar mass in galaxies 
is a crucial piece in the galaxy formation puzzle. Sub- 
stantial effort has been devoted to constraining how the 
stellar content of galaxies is related to their dark mat- 
ter halos, which provides a fundamental bridge between 
observation and theory. Such stellar mass - halo mass 
relations (M* — Mh) have yielded a variety of novel in- 
sights, including constraints on the efficiency of galaxy 
formation and connections between galaxy growth and 
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that the M* — Mh relation evolves little over most of 
cosmic history; the implications of this non-evolution are 
only beginning to be explored. 

Galaxies can be classified as 'central' or 'satellite' de- 
pending on whether they are located near the center of 
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a parent dark matter halo or whether they reside within 
a larger system. Satellite galaxies can lead tumultuous 
lives as they orbit with in the intense gravitat ional field 
of their host halo (e.g.. iKravtsov et al]|2004~bft . A num- 
ber of complex processes can therefore affect the stellar 
mass evolution of satellite galaxie s, including: (1) col d 
gas stripping due to ram pressure (IGunn fc Gotufl972lh 
(2) 'strangulation' (jLarson et al.l 1 19801 ). which removes 
the hot gas reservoir surrounding satellites, thereby re- 
ducing the fuel available to feed the cold gas disk on long 
timescales , (3) shredding of sa t ellite galaxies due to tidal 
stripping (jPurcell et al.l 120071 : IWatson et all 12012ft . and 
(4) gravitational interaction s with other nearb y galax- 
ies known as 'harassment' (|Moore et all 11998ft . These 
effects leave observational signatures on features such as 
galaxy star formation histories, colors, morphologies, and 
the amount and spatial distribution of intrahalo light. 
While much attention has focused on the evolution of 
satellites in clusters, where these effects tend to be maxi- 
mally important , the typical ~ L* satellite galaxy resides 
m poor groups flZehavi et al.ll2005t Ivan den Bosch et al] 
2007; Z ehavi et al.l 1201 lft . It may therefore be the case 
that the evolution of the typical satellite galaxy is not 
markedly different from the average central galaxy, at 
the same stellar mass. Understanding the relation be- 
tween satellite galaxies and their dark matter subhalos 



2 



(self-bound, dark matter structures orbiting in the po- 
tential of their host halo), and how this relation differs 
from the relation between central galaxies and their par- 
ent halos, has the potential to constrain the importance 
of these various complex physical processes, and is the 
primary goal of the present work. 

There are two general approaches to constraining the 
M* — Mh relation. First, there is the halo occupation 
distribution (HOD) technique (and the closely related 
conditional luminosity function technique), which adopts 
flexible functional form s for mapping gal a xies i nt o dark 
matt er halos (e.g., IFeacock fc Smith! l2000t ISeliakl 
20001: IScoccimarro et all 120011 : iBerlind fc Weinberg! 
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[Zheng et all l2007t 

Conceptually, galaxies 



are 

separated into centrals and satellites, with the centrals 
residing at rest at the center of a parent dark matter 
halo. Satellites are then modeled as orbiting within 
parent dark matter halos, with some additional assump- 
tions regarding their radial and velocity profiles. This 
approach does not explicitly use any information about 
dark matter subhalos (i.e., satellites are not constrained 
to reside within subhalos). The se cond technique 
is kn own as abundance matchi ng dKravtsov et al 
2004al IVale fc Ostrikerl 120041: [Tasitsiomi et al 
2001 IVale k Ostrikerl booa IConrov et al.1 120061: 
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In this approach galaxies 
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assigned to dark matter halos such that galaxy and halo 
mass are monotonically related with the most massive 
galaxies residing in the most massive halos. In this 
case subhalos are used to place satellite galaxies in 
the cosmic web. The abundance matching technique 
requires as an input the spatial abundance of objects, 
i.e., the luminosity or mass function of galaxies. De- 
spite its simplicity, this technique has been shown to 
reproduce a variety of observations incl uding galaxy 
two-point correlation func t ions (2PCFs) (|Conrov et all 
l2006t iMoster et all [20101: IReddick et all I20121). close 
pair counts 
2012|), M* 
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and group multiplicity functions ([Hearin et al 



20121) . However, several studies have noted discrepancies 
between abundance m atching prediction s and galaxy 
statistics at low-z fe. g..lHearin et alll2012fl and interme- 
diate redshifts (e.g., [Wetzel fc White! 120101: iGerke et all 
2012), and we expound on these shortc omings in more 
detail t hroughout the paper. Recently, iBehroozi et all 
(|2012bl hereafter B12) constrained M* - M h relations 
over most of cosmic history using a different technique, 
where satellites were placed in dark matter subhalos, 
and in which the M* — Mh relation was flexibly pa- 
rameterized and fit to a suite of data rather than being 
derived directly via the abundance matching formalism. 

A key assumption made in most modeling efforts that 
place satellites in dark matter subhalos is that the satel- 
lite galaxy - subhalo mass relation is identical to the 
central galaxy - halo mass relation (with the appropriate 
choice for the subhalo mass; see Section 2). In the present 
work we relax this assumption and attempt to separately 
constrain the satellite M* — Mh relation using galaxy 



clustering measurements and, where available, directly 
estimated satellite fr actions. This approach has also 
be taken rece ntly bvlWang et al.l (l2006h. iNeistein et all 
(|2011al) . and IRodriguez-Puebla et al.l (|2012l ). Such an 
approach can help address a variety of questions about 
satellite galaxies. For instance, How well does the stel- 
lar mass- halo mass relation of satellite galaxies trace 
the central galaxy relation? How does the satellite re- 
lation evolve with redshift? Is there strong evolution in 
the peak and characteristic mass scale for satellite galaxy 
star formation efficiency? 

The rest of this paper is organized as follows. In § [2] 
we discuss the simulation and halo catalogs used in this 
study. In |J3] we lay out our theoretical framework, in 
§[4] we present our results, and a discussion is presented 
in § [5l Finally, in § [6] w e give a summa ry of our work. 
Throughout we assume a lChabrierl (|2003l) IMF and a flat 
ACDM cosmological model with f2 m = 0.27 and fl\ — 
0.73, and H a = 70 km s" 1 Mpc -1 . 

2. SIMULATION AND HALO CATALOGS 

In the present work we employ the Bolshoi high res- 
olution collisionless N— bo dy simulation of cold dark 
matter (jKlypin et al|[201ll) . The simulation has a vol- 
ume of 250 h~ 3 Mpc 3 with 2048 3 particles and a stan- 
dard cold dark matter (ACDM) cosmological model with 
n m = 0.27, n A = 0.73, il h = 0.042, h = 0.7, a s = 0.82, 
and n s = 0.95. The simulation has a mass resolution 
of 1.9 x 10 8 Mq and force resolution of l/i _1 kpc, and 
particles were tracked from z — 80 to z — 0. Bolshoi 
was run with the Adaptive Refinement Tree Co de (ART; 
iKravtsov et~aTlll997l: IKravtsov fc Klvpinlll999D . 

Halos and subhalos were identified with the phase- 
space temporal halo finder ROCKSTAR (Be hroozi et all 
l2011alfbf ). capable of resolving halos and subhalos down 
to a maximum circular velocity of ~ 55 km s _1 . We use 
the ROCKSTAR halo catalogs at z = 0.05,0.9 and 1.9 
to model the satellite galaxy M* — Mh relation for clus- 
tering measurements at these same (median) redshifts. 
Halo masses were calculated using spherical overdensi- 
ties according to the redshift-dependent virial overden- 
sity formalism of iBrvan fc Normanl (|1998l ). 

Throughout this work we consider the halo mass def- 
inition Mp ea k for subhalos, which is the maximum mass 
that a subhalo has ever reached in its merger history. 
As we discuss in § [3j our theoretical framework is built 
off the M* — Mh relation of B12 which employs M pea k, 
thus we adopt this definition for consistency. M poa k is 
closely related to the commonly used M acc , the mass 
of a subhalo at the epoch of accretion 4 . These quan- 
tities are important because subhalos lose mass due to 
tidal stripping throughout their orbital history. How- 
ever, the galaxy is typically much more centrally concen- 
trated than the halo, and so it is expected to undergo 
tidal stripping only after its halo has been nearly de- 
stroyed. Thus we expect the galaxy content of a sub- 
halo to be more closely tied to the mass of the halo at 
the epoch of accretion or its maximum mass through- 
out its h istory, rather than the present mass of the halo 
(see e.g..lNagai fc Kravtsovll2005l: IVale fc Ostriken l2006l: 
IConrov et al.l 120061: IReddick et al.ll2012L for further dis- 

4 We have repeated our analysis using subhalos at M acc and find 
our results to be consistent with those based on M pea k- 
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Fig. 1. — "Proof of concept" demonstrating the effect of varying the satellite galaxy M* — relation on the two-point correlation 
function, satellite fraction, and stellar mass func tion. The top panel sh ows the satellite stellar mass - halo mass relation at z ~ 0, where 
the solid black curve is the Sats=Cens model of Bchroozi et al. (2012b), such that the satellite M* — relation is fixed to the central 
galaxy relation. The solid red and blue bands show the lcr deviation associated with this relation (~ 0.2 dex scatter). By populating 
a high-resolution N-body simulation with stellar masses based on this mapping, we compute the projected galaxy two-point correlation 
function for a mass threshold log(Af*) > 10.5 Mq (bottom left panel), the satellite fraction (bottom middle panel), and the galaxy stellar 
mass function (bottom right panel). Both the satellite fraction and the correlation function are very sensitive to variation in the satellite 
M* — Mft relation. In contrast, the stellar mass function is hardly affected by the ±lcr shift in the satellite M* — relation because 
satellite galaxies are a minority population. This allows us to alter the satellite relation in order to match clustering measurements and 
satellite fractions (where available) while still ensuring that the model reproduces the global stellar mass function. 
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cussion on this point). 

3. THEORETICAL FRAMEWORK 

Our approach begins with the M* — Mh relation of 
B12, which was constrained to match observational mea- 
surements of galaxy stellar mass functions (GSMFs), star 
formation rate-stellar mass relations, and the cosmic star 
formation history over the redshift range < z < 8. 
We retain the B12 relation for central galaxies, since the 
low-2; observed satellite fractions imply that one-point 
statistics like the GSMF are primarily sensitive to cen- 
tral galaxies. In our approach the satellite M* — Mh re- 
lation is constrained by observed 2PCF data and, where 
available, direct estimates of the satellite fraction from 
group catalogs. 

B12 chose a M* — Mh form that was flexible enough to 
reproduce the highly constraining low-z measurements, 
but not so flexible as to allow for over-fitting of data 
at higher redshifts where the uncertainties can be large. 
The specific functional form is given as: 

/ fM* \ \ 

log 10 (M* (M h )) = \og w (eM l ) + f (log 10 (— J J - /(0), 



(1) 



where 



f(x) = -log 10 (10 a * + l)+8 



(log 10 (1 + exp(x))F 



The virial mass at the epoch of observation is used for 
the halo mass of central galaxies, while M poa k is used for 
subhalos which harbor the satellites. The 5 parameters 
governing the fitting function are 5 :(l) the characteristic 
halo mass Mi, (2) the ratio of the characteristic stellar 
and halo masses e = Mq/M\, (3) the low-mass slope 
a, (4) the index of the subpower law 6 of the M* — Mh 
relation high-mass slope 7, and (5) the strength of the 
subpowerlaw at the massive end 8. There is also scatter 
in the M* — Mh relation, £, which is assumed to be a log- 
normal distribution and constant with (sub)halo mass. 
We will call the model where satellites and centrals are 
given the same M* — Mh relation as determined in B12, 
the Sats=Cens model. The model in which the satellite 
M* — Mh relation is allowed to differ from the central 
relation will be called the varySats model. 

The "proof of concept" of our technique is shown 
in Figure HJ The top panel shows the mean satel- 
lite galaxy M* — Mh relation at z ~ for the best-fit 
Sats=Cens model of B12, with the solid red and blue 
bands depicting the la deviation (~ 0.2 dex scatter in 
the M* —Mh relation). We populate the halos and subha- 
los from the Bolshoi simulation with stellar masses based 
on the Sats=Cens prescription and compute the pro- 
jected 2PCF (bottom left panel), the satellite fraction 
(bottom middle panel), and the GSMF (bottom right 
panel). The details of the modeling are discussed be- 
low. The crucial point is how the scatter propagates into 

5 Commonly used M* — relations have a double power- 
law form with 4 principal parameters describing a characteristic 
halo/stellar mass and low/high-mass power slopes. B12 demon- 
strated that this double power-law formalism can not accurately 
match GSMFs at various redshifts (see Appendix D of B12). 

6 A subpower law is a function that is asymptotically shallower 
than any power law, yet asymptotically greater than a logarithm. 



these three statistics. The 2PCF and the satellite frac- 
tion (/ sat = n sa t/ritot) are very sensitive to the satellite 
M*—Mh relation. For the 2PCF this is especially notable 
on scales r p < 1 Mpc (we are showing a mass thresh- 
old log(M*) > 10.5 M P) as this c orresponds to a mass 
cut used by iReddick et al.l (|2012t ) that we compare to 
throughout the paper). One can think of the 2PCF as a 
sum of two terms: on small scales, pairs of galaxies reside 
in the same host dark matter halo (the '1-halo' term), 
whereas on large scales, the individual galaxies of a pair 
reside in distinct halos (the '2-halo' term). The shape 
and amplitude of the 2PCF on small scales (r < 1 Mpc) 
are thus very sensitive to central-satellite and satellite- 
satellite pairs. / sat is sensitive to changes in the satellite 
M* — Mh relation at all galaxy stellar masses, though 
the effect is most marked at the high-mass end. How- 
ever, notice that the GSMF is hardly affected at all by 
the ±1ct shift in the satellite M*—Mh relation since satel- 
lite galaxies are subdominant by number. The fraction of 
galaxies that are satellites at low redshift reaches a maxi- 
mum vahjfirfj^ samples 
(e.g. JZheng et al. 2005; Zch avi et al.ll201lD and decreases 
with increasing redshift (pheng et al.l 120071 : iWake et al.l 
[20T1 . Therefore, we have some freedom to adjust the 
satellite M* — Mh relation to match clustering measure- 
ments while still remaining within the observational un- 
certainties of the global GSMF. 

We take the following procedure to constrain the satel- 
lite M* - M h relation: 

1. Adopt the M* — Mh relation for the central galaxies 
derived in B12. 

2. Parametrize the satellite relation such that subha- 
los are labeled according to their peak mass, M pca k, 
and pre-factors have been inserted in front of all the 
B12 M* — Mh relation parameters (as well as the 
scatter) such that, 



log 10 (Mi, sat )=Ci x log 10 (Mi) 
1 ogio(e S at)=c 2 x log 10 (e) 



= 03 x a 
- C4 x 7 
= c 5 x 8 

Csat = C 6 X £. 



Qsat 
7sat 
^sat 



(2) 



3. Start with the Sats=Cens model parameter values 
(i.e., Cj = i...6 = 1.0) and populate the halos and sub- 
halos of the Bolshoi simulation with stellar masses 
based on this M* — Mh relation. 

4. Compute the real-space two-point correlation func- 
tion from the simulation for a given stellar mass 
threshold sample. For comparison to data, convert 
either to the projected, w p (r p ), or angular correla- 
tion function, w(9). The former is computed via: 



dir. 



(3) 
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with the same upper limit of integration as used in 
the observational samples. The angular correlation 
function is computed using the Limber transforma- 
tion, 



w(9) = - / dzH{z)N(z) 2 / £(v/u 2 + D m (z) 2 6 2 )du, 
c Jo Jo 

(4) 

where c is the speed of light, H{z) is the Hubble 
constant at redshift z, N(z) is the normalized red- 
shift distribution of the galaxies in the sample, and 
D m (z) is the comoving distance to the median red- 
shift. 



5. At each redshift, compute the likelihood between 
the observed and predicted clustering for all avail- 
able stellar mass threshold samples. Incorporate 
the measured satellite fraction as an additional con- 
straint where available 7 . 

6. Use a Markov Chain Monte Carlo (MCMC) 
method, varying the parameters Ci—i-e, to com- 
pute the range of possible satellite M* — Mh re- 
lations given the data. Specifically, adopt the 
Metropolis-Hastings algorithm, which works as fol- 
lows. With the initial free parameter values from 
step 3, compute x 2 from this starting point. Steps 
are then chosen for the parameters and x 2 is com- 
puted for the new location. This new location 

or n < 



is added to the chain if \ 



2 

new 



Id 



exp[-(x ncw - Xoid)/ 2 ]' where n is a random num- 
ber between and 1. If these conditions are not 
met, then the old location is repeated in the chain. 
This process then repeats until the chain has con- 
verged 8 . 

In the present work we focus on the satellite M* — 
Mh relation at z ~ 0, 1, and 2. At z ~ we 
consider a recent measurem ent of the st e llar m ass- 
dependent 2PCF derived by iReddick et ail (|2012l) us- 
ing t he New York Univers ity Value Added Galaxy Cat- 
alog ([Blanton et al.l I2005T ). based on the Sloan Digital 
Sky Survey (SPSS) Data Rele ase 7 (jPadmanabhan et al.l 
120081 lAbazaiian et al.l I2009D 9 . We also use satellite 
fractions derived in Reddick et al. from SDSS group 
catalogs. The group catalogs were co nstructed ac- 
cordin g to the methodology described in I Tinker et al.l 
(|2013l) . At z ~ 1 we employ two-point clustering mea- 
surements from the DE EP2 Galaxy Redshi ft Survey 
(jNewman et al.ll2012D by iMostek et al.l (|2012l , hereafter 
M12). At z ~ 2 constraints are provided by angular 
clustering measurements from the NEWFIRM Medium 

7 Direct estimates of satellite fractions are not available at z ~ 1 
and 2 and so this additional constraint is not considered at those 
epochs. It is also true that w p (r p ) and / sa t are highly covariant. 
We do not take this into account in our modeling, though this 
should be incorporated in future work. 

8 We incorporate the full covariance matrices for our model fit- 
ting at z ~ and 2. Covariance matrices are not available for the 
DEEP2 z ~ 1 samples. We also assume no priors on the parameters 
and we start numerous chains at a wide range of initial c pre-factor 
values to confirm that they converge to the same parameter values . 
For more details on MCMC techniques, see Dunklcv et al. (2005). 

9 We have converted all of the Reddick et al. (2012) results from 
H = 100 h km s" 1 Mpc" 1 to H = 70 km s _1 Mpc -1 . 



Band S urvey (NMBS; Ivan Dokkum et alj |2009t) as de- 
rived in I Wake et al.l In the following section we 
confront our model with these data. 

4. RESULTS 

We now compare the Sats=Cens and varySats models 
to observational constraints at z ~ and 1, and the 
Sats=Cens model to data at z ~ 2. We derive satellite 
M* — Mh relations at < z < 1 and investigate its 
evolution and how it compares to that of central galaxies. 
We also present satellite fractions over the interval < 
z < 2. 

4.1. Clustering and the Satellite M* — Mh Relation at 

z ~ 

Figure [2] shows the projected 2PCF at z ~ at 
three separate stellar mass thresholds: log(M*) > 10.1, 
10.5, and 11.1 Mq. The black dashed curves show the 
Sats=Cens model where the clustering is computed from 
the best-fit M* — Mh relation given by B12. We empha- 
size that this curve is not a fit to the clustering data, 
rather it is a direct prediction of the B12 M* — Mh re- 
lation. The Sats—Cens model provides an excellent fit 
to the data at the largest stellar mass thresholds and on 
large scales, but under-predicts the clustering on small 
scales in the two lower mass samples. 

This under-prediction is probably not due to reso- 
lution effects in t he Bolshoi simulation , as discussed 
in Appendix B of Reddick 'et al.l (120121). In ad dition, 



IWatson et all (|2012l ) used the lZentner et al.1 (|2005| ) semi 
analytic model for subhalos, which has the capabil- 
ity of tracking subhalos down to circular velocities ~ 
km s _1 , to test resolution effects in numerical simu- 
lations. They performed a test in which they selected 
a > 210 km s _1 threshold (corresponding to ~ L* 
galaxies and brighter) with two additional v max thresh- 
olds of v max > 20 km s" 1 and v max > 80 km s _1 to 
mimic the effect of a resolution limit of a given simu- 
lation. Increasing this threshold from 20 to 80 km s _1 
resulted in a maximum ~ 20% decrease in the 2PCF 
at scales r < IMpc. Increasing from 20 to 40 km s _1 
only had a few percent effect on the small-scale ampli- 
tude. Though not explicitly shown in their paper, they 
repeated this test for v^ ax > 170 km s _1 and found 
nearly the same result, implying that the resolution limit 
of 55 km s _1 should not have a significant effect on the 
2PCF for the lowest stellar mass threshold considered in 
Figure [21 

We have also constructed 1000 Sats—Cens models 
drawing from the errors in the B12 best- fit parameters 
and found the la distribution in the 1000 2PCFs to still 
be significantly discrepant from the data on small scales 
at lower masses. Thus, no Sats=Cens model can provide 
a satisfactory fit to the small-scale clustering for the low- 
est stellar mass threshold sample. 

The Sats=Cens model succeeds on large scales because 
the large-scale clustering is only weakly sensitive to the 
treatment of satellite galaxies. On smaller scales the 
2PCF is much more sensitive to the relation between 
satellite galaxy mass and halo mass. The varySats model 
allows for freedom in the satellite M* — Mh relation and 
thus is able to more closely match the data shown in Fig- 
ure [U though there is clearly still unresolved tension at 
the 1- and 2-halo junction at ~ 1 Mpc. 
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0.1 1.0 10.0 0.1 1.0 10.0 0.1 1.0 10.0 

r [Mpc] r [Mpc] r [Mpc] 

Fig. 2. — Comparison between model a nd ob served projected correlation function as a function of stellar mass for three SDSS threshold 
samples as measured by Roddick et al. (2012): log(Af*) > 10.1, 10.5, 11.1 Mq . The black dashed curve shows the Sats=Cens model 
prediction where the satellite M* — relation is fixed to that of the central galaxies. The red solid curve shows our best-fit model in 
which the M* — relation for satellite galaxies is allowed to vary to better match observed clustering and satellite fractions (see Figurc[3]L 
There is still tension between the varySats model and the data at the 1-2 halo transition, especially at lower masses. 
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Fig. 3. — Top panel: Best-fit model satellite fraction (red curve) 
as a function of galax y stellar ma ss (black data points with error 
bars are from Rcddick ct al. 2012). The Sats=Cens model clearly 
fails to reproduce the measured satellite fractions at low masses, 
while the varySats model is in good agreement with the data. Bot- 
tom panel: Residuals between th e observed an d mod el galaxy stel- 
lar mass function; data are from lBaldrv et al.l 12008T ). The galaxy 
stellar mass function is dominated by central galaxies, so varying 
the satellite galaxy M * — M h relation has little effect on the galaxy 
stellar mass function and is within the errors. 



In addition to the 2PCF, the varySats model is also 
constrained by the measured satellite fraction, shown in 
Figure [3] As we see in the top panel of Figure [31 the data 
require a larger satellite fraction for log(M*) < 11.0 M Q 
than the prediction of the Sats=Cens model. The bot- 
tom panel illustrates that the changes to the global 



GSMF of IBaldrv et all (|200l (the same GSMF used in 
B12) implied by the varySats model are small and within 
the error bars. 

To evaluate whether (and to what degree) the varySats 
model is preferred over the Sats = Cens model we em- 
ploy the Aikake Information Criterion (AIC; see Liddle 
120071 for details). AIC is preferred over a reduced x 2 
comparison since we do not take into account the co- 
variance between / sat and w p (r p ), which makes it diffi- 
cult to discuss the goodness-of-fit in an absolute sense. 
Instead, using AIC we are able to comment on model 
preference. AIC is defined as AIC= 2k + Xmin> where 
k is the number of free parameters in the model and 
in our case x 2 = X 2 (/sat) + X 2 ( w p( r p))- It is stressed 
in the literature to add a correction factor to AIC to 
take into account small sample size, such that AIC C = 
AIC + 2k(k + l)/(N - k - 1), where N is the number 
of data points (N — 44 for the z ~ case). For the 
Sats — Cens model, k = 0, X% a ts=Cena = 220.1, the cor- 
rection factor is 0, so AIC C c . „ — 220.1. For the 

' L - bats— Cens 

varySats model, k — 6, xlarySats = 68.5, the correction 
factor is 2.2, thus AIC C c , = 82.7. When judged on 

1 "-'vary bats J ° 

Jeffreys scale for interpretation, the AAIC C implies that 
the varySats model preference is considered 'decisive'. 

Turning to Figure El the top panel shows the satellite 
M* — Mh relation: satellite galaxy stellar mass vs. sub- 
halo peak mass, M pca k. Again, the dashed black curve is 
the Sats=Cens model of B12, and the red curve shows 
the best-fit result from the varySats model. The solid red 
band shows the la uncertainty in our varySats model. 
The varySats model implies that above the complete- 
ness limit of log(M*) > 10.1 M Q of lReddick et all ()2012D 
there should be, on average, more stellar mass in satellite 
galaxies than the Sats=Cens model prediction by ~ 10% 
at all subhalo masses where the model can be constrained 
(although the two relations are consistent at the 2 — 3cr 
level for subhalo masses > 10 13 /i -1 Mq). This implies 
that there will be more satellite galaxies in each SDSS 
stellar mass w p (r p ) threshold sample, which will boost 
the clustering and the satellite fractions. 
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Fig. 4. — The satellite galaxy M* — M h relation at z ~ 0. Top 
panel: The black curve represents the Sats=Cens model. The thin 
solid lines demarcate the mass range over which the model has 
been constrained by data. The red curve shows the best-fit rela- 
tion for the varySats model. The lighter red band represents the \a 
uncertainty in the model. The varySats model is driven to higher 
masses in order to boost the clustering and satellite fractions. Bot- 
tom panel: Ratio between stellar mass and halo mass as a function 
of halo mass (the integrated star formation efficiency). The peak 
in star formation efficiency and the characteristic (sub)halo mass 
is nearly identical in both cases highlighting the mild stellar mass 
evolution in satellite galaxies and their similarity to the known lack 
of evolution of central galaxies. 



The bottom panel shows the fraction of available 
baryons that have turned into stars (the integrated star 
formation efficiency) to illustrate at what subhalo mass 
galaxy formation is most efficient. Above the stellar mass 
limit where the model can be constrained, star formation 
in the varySats model has been slightly more efficient, 
on average, at all subhalo masses than central galaxies 
of the same mass. On the other hand, the characteristic 
mass for peak satellite galaxy star formation is approxi- 
mately the same for central and satellite galaxies, though 
we urge caution with this conclusion as we are unable to 
probe the full range of the "turn over" in the M* — Mh 
relation. 



4.2. Clustering and the Satellite M* 

z ~ 1 



Mh Relation at 



As we turn to the z ~ 1 results, we find that the 
Sats=Cens model is in rather good agreement with the 
z ~ 1 DEEP2 data as shown in Figure [5] Again, 
the black dashed curve shows the Sats=Cens model 



prediction and the red curve is our best-fit varySats 
model with the light red band representing the la un- 
certainty in the model. Satellite fraction measurements 
were not available fo r these data, however, we use the 
iPerez-Gonzalez et al.l (|2008h total GSMF to put bounds 
on allowable satellite M* — Mh relations. Again, this is 
the same GSMF data used in B12. Our varySats model 
essentially reproduces the Sats=Cens model for the clus- 
tering at log(M*) > 10.5 M Q and is slightly boosted 
at all physical separation scales of the data for the 
log(M*) > 10.8 Mq sample. As was done in § 14.1 1 we can 
compute AIC C to compare the performance of the two 
models. In this case, x 2 = X 2 (GSMF) + X 2 (%>(?>)) and 
N = 45, with X 2 sat s =cens = U8-2 and xLrySats = 93.9. 
The straightforward AAIC C computation indicates that 
there is no decisive preference for the varySats model 
over the Sats=Cens model. 

In Figure |6] we show the best-fit M* — Mh relations. 
Here the errors on the satellite relation are larger than 
at z ~ 0, both because the clustering data carry larger 
errors and because we do not include measured satellite 
fractions in the fitting. Nonetheless, it is clear that the 
satellite relation is very similar to the central relation 
at z ~ 1. Future data from larger area surveys at z ~ 1 
should provide stronger constraints on the satellite M* — 
M h . 

4.3. Clustering at z ~ 2 

IWake et"ail (|2011h recently presented stellar-mass de- 
pendent angular correlation functions of galaxies at z ~ 2 
based on the NMBS. The errors on the measurements 
are relatively large (compared to lower redshift) owing 
to the relatively small area of the NMBS fields. For this 
reason we did not attempt to fit the VarySats model to 
the data. Instead we simply compared the Sats=Cens 
model prediction to the existing clustering data. The 
result is shown in Figure [7] Again we emphasize that 
the Sats=Cens model makes a prediction for the clus- 
tering of galaxies, as this model was only tuned to fit 
the z ~ 2 GSMF. The model fits the data remarkably 
well over all scales and all mass ranges probed by the 
data. Future surveys covering a larger volume will be 
necessary to make more quantitative statements regard- 
ing the similarity of the satellite and central M* — Mh 
relations. 

4.4. Satellite Fractions at z ~ 0, 1 and 2 

Our modeling describes how galaxy stellar mass should 
be mapped to host halos and subhalos. This mapping 
can be used to infer satellite fractions as a function of 
stellar mass and redshift. Figure [5] shows the satellite 
fraction where constraining data exists as a function of 
galaxy stellar mass from our best-fit varySats model at 
z ~ and I (red and blue curves). The green curve is 
the Sats=Cens model at z ~ 2. The solid bands are 
the associated la errors. The varySats model was not 
considered at z ~ 2, thus no errors are shown for this 
curve. Satellite fractions decline wit h increasing redshift 



in agreement with previous studies ( Zentner et al.l 120051: 
Zheng et aJJ l2fK)7t iWetzel fc Whit^ 1201(1 iWatson et al.1 
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Fig. 6. — Same as Figure [4] but for z ~ 1. The errors on the 
Z ~ 1 varySats model are larger than at z ~ 0, with the result that 
the M* — Mh relation for this model is consistent with the relation 
for the simple Sats=Cens model. 



5.1. Implications for Satellite Galaxy Evolution 

The primary result of this paper is that the mean re- 
lation between galaxy stellar mass and halo mass is very 
similar for central and satellite galaxies since at least 
z ~ 2. At z ~ there is a slight difference in the sense 
that satellites appear to have somewhat higher stellar 
masses at fixed halo mass, while at z ~ 1 — 2 there is 
no evidence for a difference between the satellites and 
central galaxy M* — Mh relation. In this section we ex- 
plore the implications of these results for the evolution 
of satellite galaxies. In order to focus the discussion, in 
Figure [5] we combine into a single diagram the central 
and satellite relations over the interval < z < 2. The 
central galaxy relations (dashed curves) are the global re- 
lations from B12 and the satellite relations (solid curves) 
are those derived herein. 

There is an emerging consensus that the global 

M* — Mh relation evolves little if at all sinc e z = 

2 (see Figure M and also e.g. , lYang et ail 120111; 
Moster et all l2012t iBehroozi et al.l l2012bt IWang et al.l 
2012t iLeauthaud et al.l l2012f ) . As discussed in previous 
sections, the satellite fraction in mass-limited samples 
is always < 30% and decreases with increasing redshift. 
The global relation is therefore, to a good approxima- 
tion, the central galaxy relation. The weakly evolving 
M* — Mh relation for central galaxies is one of the main 
facts that explains the lack of evolution in the satellite 
galaxy M* — Mh relation. In the discussion that follows 
we will assume for simplicity that the central relation is 
constant in time. 

First, it is worth briefly considering scenarios that are 
inconsistent with our results. Satellite galaxies, on aver- 
age, cannot have experienced large (order unity) growth 
in stellar mass since their accretion, nor can they have ex- 
perienced a large loss in mass (unless they are altogether 
destroyed). These possibilities would result in satellite 
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M * — Mh relations significantly offset from the central 
galaxy relation. Of course, since the resulting satellite 
relation is an integral over the past history of the satel- 
lite, in principle our results allow for a scenario in which 
the average satellite has both gained substantial mass via 
e.g., star formation and lost substantial mass via e.g., 
tidal stripping. This seems implausible, at least for the 
average satellite. It is therefore clear that any successful 
scenario must allow for only modest growth in the stellar 
mass of satellite galaxies while they are satellites. 

In another scenario for satellite galaxy evolution, stel- 
lar mass growth is immediately halted when a galaxy 
becomes a satellite. In this case the satellite M * — Mh 
relation would then be identical to the unevolving cen- 
tral relation. This would be so because for any scenario 
the subhalo peak mass, M pea k, does not evolve after ac- 
cretion by definition, and in this scenario neither does 
M*. Such a possibility is reminiscent of previous gener- 
ations of semi-analytic models, in which satellite SFRs 



were truncated shortly after infall. However, these ear- 
lier models have been shown to overproduce the frac- 
tion of quiescent sate l lites as a functio n of mass (e.g., 
IWeinmann et alj|2006t front et al.ll2008t) . 

We now turn to a more plausible scenario for satellite 
galaxy evolution. In this scenario, all satellites continue 
to grow in stellar mass as if they were still central galax- 
ies for a significant fraction of their lifetime (see also 
iWetzel et al.ll2012bl ). If we adopt the relation between 
star formation rate (SFR) and s t ellar mass measured in 
the local universe bv lSalim"e t al. (2007|) along with a typ- 
ical time since accretion for satellites of 4 Gyr (see e.g., 
iZentner et all 120051; IWetzel et alj|2012b| ). then the aver- 
age galaxy with M* ~ 10 11 M will have increased its 
mass by ~10% since becoming a satellite, if it evolved 
like a central galaxy for its entire lifetime. Owing to the 
fact that the SFR/M* relation is so shallow, the frac- 
tion of mass grown is a weak function of mass. Thus, in 
this scenario, the satellite M* — relation would be 
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in bins of Alog(M*) = 0.1 as inferred from our best-fit varySats 
model at z ~ and 1 and from the Sats=Cens model at z ~ 2. 
Satellite fractions are only plotted over the mass range where there 
is constraining data. 



shifted approximately 10% higher compared to the cen- 
tral relation (i.e., at fixed halo mass satellites would have 
~ 10% more stellar mass). It is interesting that this shift 
is similar to what we find herein at z ~ 0. Given the sys- 
tematic uncertainties that exist at the 10% level in these 
relations we caution against over-interpreting this simi- 
larity, although we will now briefly explore why satellites 
may evolve like centrals for most of their evolution. 

Satellite galaxies can be subject to numerous physi- 
cal processes capable of altering their evolution. This 
has been well-documented in very dense environments, 
such as large galaxy groups and clusters. For instance, 
ram pressure stripping and strangulation can reduce the 
gas supply, stifling star formation. However, these pro- 
cesses are less effective in modest-sized groups where 
most satellites reside. In a given mass-limited sample of 
galaxies, most satellites are found in relatively quiescent 
environments, e.g., in a halo that contains one central 
and one satellite gala xy. This is evident in the group 
multi plicity function (|Berlind et al.l 120061: iHearin et all 
20121) and in halo occupat i on models (IZehavi et alj|2005t 
van den Bosch et alJl2007t IZehavi et al.ll2011f ). Thus, for 
the typical satellite, physical processes capable of alter- 
ing its evolution are relatively inefficient. 

Perhaps more importantly, dynamical friction 
timescales at the group scale are shorter than at the 
cluster scale (for a satellite of the same mass), meaning 
that satellite galaxies in groups spiral toward the center 
and merge with the central galaxy more quickly than in 
rich clusters. This implies that there is a relatively small 
amount of time for the group environment to shape the 
physical properties of a satellite before that satellite is 
destroyed and/or merges with the central object. This 
effect becomes even stronger at higher redshifts because 
dynamical times are shorter. This also partly explains 
the lower satellite fractions at higher redshifts discussed 
in Section IP1 

This picture was motivated in part by the recent re- 
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relations at z ~ 0, 1 and 2. 



suits of iWetzel et al.l (|2012bl ). who analyzed the proper- 
ties of low redshift satellite and central galaxies including 
stellar masses and star formation rates. These authors 
concluded that, after it is accreted, the typical satellite 
continues to form stars as if it were still a central galaxy 
for 2—4 Gyr before star formation is quenched. This con- 
clusion follows naturally from the fact that there exists 
a strong bimodality in satellite galaxy SFRs combined 
with the fact that the locations of the two peaks are sim- 
ilar to those of centra l galaxies with the same stellar mass 
(jWetzel et al.ll2012aT ). 

The similarity between satellite and central M * — Mh 
relations thus appears to be a natural outcome of sev- 
eral facts: 1) the weakly or unevolving M* — Mh relation 
for central galaxies; 2) the hierarchical growth of struc- 
ture, which ensures that the majority of satellite galaxies 
(at least with M* > 10 10 M o ) are not strongly affected 
by their environments for most of their relatively short 
lives; 3) typical SFRs at late times are low, resulting 
in long mass doubling times; and 4) the combination of 
the relatively short satellite lifetimes and long mass dou- 
bling times suggests that the typical satellite will not 
have grown substantially in mass during its lifetime, and 
therefore the satellite relation will not differ significantly 
from the central galaxy relation when the satellite was 
accreted, which in turn differs little from the present day 
central relation. 

5.2. Comparison to Previous Work 

There have been several previous studies that at- 
tempted to separately constrai n the satelli t e and cen- 
tral galaxy M* - M h relations. IWang et all (|2006l ) were 
the first to consider such a distinction. They adopted 
a method that used galaxy positions and velocities de- 
termined from the orbital and merging histories of ha- 
los and subhalos from a full cosmological simulation (as 
opposed to a statistical HOD approach) to predict ob- 
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served galaxy statistics. Galaxy properties such as lu- 
minosity and stellar mass were parameterized as a func- 
tion of M acc , where the functional forms of the L — Mh 
and M* — Mh relations were adopted from a semi- 
analytic model as reference. They were able to success- 
fully model the observed stellar mass/luminosity func- 
tions and 2PCFs at low redshift and found central and 
satellite M* — Mh relations that were identical within 
the errors. These results are broadly consistent with 
our work, although recent high-precision measurements 
of w p (r p ) and / Ba t argue for a slight difference between 
t he central and satellite M* — M h relations. 

iRodriguez-Puebla et a l. (2012) also considered distinct 
M* — Mh relations for centrals and satellites in an 
abundance matching context. They used the z ~ 
GSMF decomposed into central and satellite galaxies and 
used abundance matching separately for the centrals and 
satellites. They found models based on M acc were in 
agreement with measured conditional stellar mass func- 
tions and correlation functions only when the satellite 
M* — Mh relation was slightly different from the central 
relation (by < 0.1 dex at M > 1O 12 M ), a result in 
g eneral accord with our low-z findings. 

iNeistein fc Khochfarl (|2012f ) used abundance match- 
ing within a full cosmological simulation to further 
investigate th e M* — Mh connection. Motivated by 
INeistein et al.1 ()2011al ). satellite stellar mass depended 
on both M acc and the host halo mass at z = 0, and the 
satellite M* — Mh relation was allowed to differ from the 
central relation. The authors considered the abundance 
of galaxies, the 2PCF and weak lensing measurements in 
their modeling and they found that the M* —Mh relation 
for satellites at z — was poorly constrained, such that 
a wide range of relations were consistent with the data, 
resulting in satellite fractions spanning / sat = 0.05 — 0.8. 
This result suggested that the inclusion of the measured 
satellite fraction as an additional statistic could be key 
to further constraining the satellite M* — Mh relation. 

Each of these studies attempted to separately constrain 
the satellite and central galaxy M* — Mh relations at low 
redshift, allowing us to piece together a clearer picture 
of the two relations. While there are several differences 
in the methodology between our work and the aforemen- 
tioned studies, the principal distin ction is our use of the 
measured satellite fractions from iReddick et al.l ()2012f ) 
as an additional model constraint. The power of using 
satellite fractions to constrain the satellite M* — Mh re- 
lation can also be seen in our Figure 1. However, at 
z ~ 1 satellite fractions are not available, and our de- 
rived satellite M* — Mh relation is consequently more 
uncertain. In fact, the uncertainties we derive are qual- 
itatively consistent with th e error envelope quoted by 
INeistein fc Khochfarl (|2012t ) based on their analysis of 
z ~ data, also without satellite f ractions. Addition- 
ally, one can see from Figure 3 of IWang et al.l (120121) 
that t he successful modeling of 2PCFs from lWang et al.l 
(2006) ove r-predicts sa t ellite fractions from the group 
catalogs of lYang et al.l (|2008l ). Our results argue that 
the observed 2PCF and GSMF are not enough to pro- 
vide strong constraints on the satellite M* — Mh relation. 
However, estimates of the mass-dependent satellite frac- 
tion, in conjunction with these other statistics, result in 
much smaller errors on the inferred satellite M* — Mh 
relation. 



It is also important to emphasize here that we do not 
allow freedom in the central relation in order to limit 
computation time. We argue that this approximation 
is justified because the global GSMF and 2PCFs on 
large scales (where central galaxies dominate) are well-fit 
by our model. However, incorporating a parameterized 
M* — Mh central relation may result in a larger range of 
possible models for the satellites. We leave this possibil- 
ity for fut ure work. 

Finally, IReddick et al.l (|2012j ) showed that v pea k, the 
maximum circular velocity analog to M pca k, was the best 
halo property for matching low-z clustering data. We 
have adopted the latter herein in order to exploit the 
newly calibrated M * — Mh from B12, which was derived 
for masses rather than circular velocities. It remains an 
open question if the satellite and central relations would 
be even more similar if circular velocities were used in 
place of halo masses. 

6. SUMMARY 

In this paper we have explored the relation between 
satellite galaxies and their dark matter halos since z ~ 2. 
Global M* — M h relations have been estimated out 
to high redshi f t (lYang et all l2011t iMoster et al.1 120121 : 
iBehroozi et"aH l2012bt IWang et al.l (2012J) by adopting 
functional forms for how galaxies populate dark mat- 
ter halos in a manner flexible enough to simultaneously 
match GSMF measurements at multiple redshifts. While 
this can yield valuable insight into galaxy formation pro- 
cesses, it is principally constraining the evolution of stel- 
lar mass in central galaxies, since satellite galaxies are a 
minor contribution to GSMFs. Incorporating recent pre- 
cision measurements of the galaxy 2PCF as a function of 
stellar mass at z ~ 0, 1 and 2, a statistic which is highly 
sensitive to the satellite galaxy population, introduces an 
additional constraint to better understand the evolution 
of stellar mass in satellite galaxies. To model the satellite 
galaxies, we started by adopting the M* — Mh relation 
of B12 which was shown to accurately reproduce GSMFs 
out to z ~ 8. We applied this relation to halos and sub- 
halos in the Bolshoi dissipationless N— body simulation 
identified by the ROCKSTAR halo finder at z ~ 0, 1 and 
2. We then calculated the projected 2PCFs to compare 
to observations at z ~ and 1 as well as angular corre- 
lation functions at z ~ 2. This approach was named the 
Sats=Cens model as the B12 M* — Mh relation gives the 
same functional form for both central and satellite galax- 
ies. The Sats=Cens model was shown to slightly under- 
predict satellite fractions and clustering measurements 
in low-mass galaxies at low redshift. This prompted us 
to introduce a model that allowed for more flexibility in 
the satellite M* — Mh relation, the varySats model. We 
retained the B12 best- fit relation for centrals and fit for 
the satellite relation, using an MCMC method to probe 
parameter space. Our technique is essentially a hybrid 
between the halo occupation distribution and abundance 
matching methods, in the sense that we adopt the abun- 
dance matching requirement that satellites reside within 
subhalos while allowing for a flexibility in the galaxy- 
halo relation common to halo occupation models. Our 
modeling approach yields the following principal results. 

• The clustering and satellite fraction data at z ~ 
imply that the satellite galaxy M* — Mh relation is 



12 



nearly identical to the central galaxy relation, with 
a preference for the stellar mass of satellites to be ~ 
10% larger than centrals at fixed halo mass (the two 
relations are consistent at 2 — 3<r for subhalo masses 
> 10 13 M Q ). At z ~ 1 the satellite and central 
M* — Mh relations are indistinguishable within the 
errors. 

• Our models provide an excellent fit to the stellar- 
mass dependent clustering of galaxies at z ~ 
1. This stands in contrast to several previous 
studies that were unable to reproduce the ob- 
served z ~ 1 £?-band luminosity-d ependent cluster- 
ing through abundance matching (|Wetzel fc White! 
20101 IGerke etalllMtl . The success herein may 
perhaps be due to the fact that stellar mass is more 
strongly correlated with dark matter halo mass 
than £?-band luminosity. 

• At z ~ 2, we find that a simple model in which the 
satellite M* — Mh relation is the same as the cen- 
tral relation is in good agreement with two-point 
clustering data at all scales and all stellar masses. 
However, at present the errors on the data are too 
large to constrain the more flexible satellite model. 

While the detailed physics of galaxy formation is com- 
plex and highly nonlinear, our results add to the emerg- 
ing consensus that the basic properties of galaxies are 
surprisingly regular and simply connected to the un- 



derlying dark matter distribution. The M* — Mh rela- 
tion for central g alaxies has been s hown to evolve lit- 
tle since z ~ 4 (lYang et all 120111; iMoster et ail 120121; 



IBehroozi et al.l l2012b ). Alo ng with this lack of evolu- 
tion. iBehroozi et al.l ( 2012al ) showed that there exists a 
nearly constant peak in star formation efficiency in halos 
with characteristic mass 10 11 7 since z ~ 4, consistent 
with theoretical predictions of t he halo mass scale for 
the shock heat ing of accreted gas dRees fc Ostrikerlll977|; 
White & R"eeill978t iBlumenthal et al.lll984l;lKeres et ail 
20051; iDekel fc Birnboimll2006t) . These results for central 
galaxies also appear to hold for the satellites. The reason 
for this is not yet entirely clear, but it is likely related to 
the fact that the typical ~ L* satellite galaxy resides in 
a poor group where transformation processes are weak 
and lifetimes are short. 
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